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The cosmic microwave background (CMB) energy spectrum is a near-perfect blackbody. The 
standard model of cosmology predicts small spectral distortions to this form, but no such distortion 
of the sky-averaged CMB spectrum has yet been measured. We calculate the largest expected 
distortion, which arises from the inverse Compton scattering of CMB photons off hot, free electrons, 
known as the thermal Sunyaev-Zel’dovich (tSZ) effect. We show that the predicted signal is roughly 
one order of magnitude below the current bound from the COBE-FIRAS experiment, but can 
be detected at enormous significance (> 1000a) by the proposed Primordial Inflation Explorer 
{PIXIE). Although cosmic variance reduces the effective signal-to-noise to 230a, this measurement 
will still yield a sub-percent constraint on the total thermal energy of electrons in the observable 
universe. Furthermore, we show that PIXIE can detect subtle relativistic effects in the sky-averaged 
tSZ signal at 30 (t, which directly probe moments of the optical depth-weighted intracluster medium 
electron temperature distribution. These effects break the degeneracy between the electron density 
and temperature in the mean tSZ signal, allowing a direct inference of the mean baryon density 
at low redshift. Future spectral distortion probes will thus determine the global thermodynamic 
properties of ionized gas in the universe with unprecedented precision. These measurements will 
impose a fundamental “integral constraint” on models of galaxy formation and the injection of 
feedback energy over cosmic time. 

PACS numbers: 98.80.-k, 98.70.Vc 


Introduction — The energy spectrum of the cosmic mi¬ 
crowave background (CMB) radiation is extremely close 
to a perfect blackbody mm- However, distortions to 
this blackbody spectrum arise from physical processes 
in both the early (redshift z > 1100) and late {z < 
1100) universe. These processes include energy injec¬ 
tion from decaying or annihilating particles (e.g., mm), 
the dissipation of small-scale primordial density fluctua¬ 
tions (e.g., [5H7]), and more exotic processes, such as en¬ 
ergy injection from cosmic strings [8] or primordial black 
holes [S] . The most well-understood distortion is that due 
to the thermal Sunyaev-Zel’dovich (tSZ) effect, which is 
the inverse Compton scattering of CMB photons off hot, 
free electrons [101 El] ■ This signal is predominantly due 
to hot gas in the intracluster medium (ICM) of galaxy 
groups and clusters at z <2. 

Many analyses have constrained cosmology and ICM 
physics using tSZ measurements of cluster samples or 
the fluctuation properties of the tSZ field (e.g., (Mni). 
However, the mean tSZ signal of the universe has yet 
to be detected — in fact, no global spectral distortion 
of the CMB has been seen to date. The tightest con¬ 
straint comes from the COBE-FIRAS experiment, which 
found |(y)| < 1.5 x 10“^ at 95% confidence [2|, where 
y is the Compton-y parameter characterizing the (non- 
relativistic) tSZ distortion, 

y(n) =[ dlne{n,l)kBTe{h,l). (1) 

meC^ J 

Here, and Tg are the electron number density and tem¬ 


perature, respectively, ctt is the Thomson cross-section, 
mgC^ is the electron rest-mass energy, and fee is Boltz¬ 
mann’s constant. This expression assumes Tg 3> Tcmb, 
which is valid for all sources in our analysis. Here, 
JcMB = 2.726±0.001 K is the average CMB temperature 
today [IH]. The integral is taken along the line-of-sight 
distance I, while n is a position on the sky. Angle brack¬ 
ets (• • •) refer to averages taken over the sky. We show 
below that the expected (y) signal is roughly one order 
of magnitude below the COBE-FIRAS limit and will be 
easily detectable with currently proposed experiments. 

Recently, interest has grown in new measurements of 
the CMB energy spectrum [BEn]. Such measurements 
have the potential to detect or place interesting con¬ 
straints on a range of new physics, such as primordial 
non-Gaussianity m or decaying dark matter [B 1^ - 
The CMB specific intensity, at frequency iz in a 

direction on the sky, h, can be decomposed into: 

= H,(TcMB + Ar) + AJj(n) + A/f2(h) + 
AU(h) + A/°*h“(n), (2) 

where Bt,{T) = ^^/(e“ — 1) with x = hv/^k-aT). 
The focus of this paper is the sky-averaged tSZ signal, 
(A/*®^). Since the exact value of the primordial CMB 
temperature is not known to the nK level accessible to 
future experiments, it must be determined in the analy¬ 
sis, which is represented by the AT term in Eq. ([^. We 
assume a fiducial value of AT = 1.2 x 10“"^ K, within 
the Ict error from COBE-FIRAS. The A/J (n) term in 
Eq. (1) represents anisotropies in the primordial CMB 
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temperature, which in a given pixel are indistinguishable 
from AT. The other terms in Eq. ([^ represent contri¬ 
butions from the /i-distortion (A/^) and all other con¬ 
tributions to the sky intensity such as other 

sources of CMB spectral distortions or foregrounds (e.g., 
dust or synchrotron emission). For the ^-distortion, we 
assume a constant amplitude of /r = 2 x 10“®, consis¬ 
tent with predictions from the damping of small-scale 
modes in the concordance ACDM cosmology Pll ■ We set 
A/°‘*'er _ Under these assumptions, we examine de¬ 
tection prospects for (A/*®^) for upcoming experiments, 
specifically the proposed Primordial Inflation Explorer 
{PIXIE) [T9]. 

We show that PIXIE can make a high-significance 
measurement of not only (y), but also the mean opti¬ 
cal depth-weighted ICM electron temperature, {k^Tflir^ 
where T(n) = ctt / dlnflnfl). The {k-BTe)r sensitivity 
arises from the relativistic tSZ signal due to the hot ICM 
of galaxy groups and clusters, and thus calibrates the 
amplitude of the cluster mass - temperature scaling rela¬ 
tion that will be fundamental for cosmological constraints 
from X-ray cluster counts in the upcoming eROSITA mis¬ 
sion |5S|. Measurements of (y) and {kBTe)r provide “inte¬ 
gral constraints” on models of galaxy formation over cos¬ 
mic time, similar to the constraint on the optical depth to 
reionization, tcmBi provided by large-scale CMB polar¬ 
ization measurements. CMB spectral distortion signals 
thus constrain the uncertain but crucial feedback mech¬ 
anisms used in hydrodynamical simulations of cosmic 
structure formation (e.g., [2M28] b Moreover, constraints 
on the electron temperature through the relativistic tSZ 
effects allow an inference of the mean electron density 
from (y), yielding a definitive answer to the long-standing 
“missing baryons” problem at low redshift (e.g., 

In a companion paper |32j (hereafter B15), we provide 
technical details of our model and calculations, and com¬ 
bine the forecasted measurements of (y) and {k^Tflir with 
other data to constrain cosmological parameters and ICM 
gas physics models. 

Theory — The sky-averaged tSZ signal receives contri¬ 
butions from three main components: the ICM of col¬ 
lapsed halos (galaxies, groups, and clusters), the inter- 
galactic medium (IGM) between halos, and the epoch of 
reionization. The latter two contributions are subdomi¬ 
nant to that from the ICM, but all are included in our 
analysis. Physically, the global tSZ signal probes the in¬ 
jection of energy into ionized gas over the history of the 
universe. Using Eq. 0 . 

where the integrals are taken over the observable uni¬ 
verse, Pe = n^k^Tf, is the electron pressure, and = 

^coii_^cooi_|_^inj jg total thermal energy in electrons. 
Here, E™" is the energy in electrons due to gravitational 
collapse during structure formation, is the electron 


energy lost to cooling processes, and E™^ is the energy 
injected into electrons by feedback processes (e.g., from 
supernovae and active galactic nuclei). For sufficiently 
high electron temperatures (feeTe/^eC^ ^ 10“^), rela¬ 
tivistic corrections to the tSZ effect become important, 
inducing an additional dependence on Te [55H55] . 

The ICM and IGM contributions to (AJ*®^) can be 
calculated directly from hydrodynamical simulations [36l - 
[5^ . The ICM signal can also be calculated in the halo 
model (e.g., [ID]) using calibrated prescriptions for the 
electron pressure and temperature as a function of halo 
mass and redshift (Ml HD ■ We define the IGM as gas 
located at greater than 2.5 times the virial radius of any 
collapsed halo. 

The sky-averaged tSZ signal from ICM electrons is 
(e.g., Eq. (2) of Ref. |12|) 



(4) 


where d^V/dzdH is the comoving volume per steradian 
at redshift z, dn/dM is the halo mass function (the co¬ 
moving number density of halos as a function of mass M 
and redshift), and A/*®^(n, M, z) is the tSZ signal at an¬ 
gular position n on the sky with respect to the center of 
a cluster of mass M at redshift z. In the non-relativistic 
approximation, 

A/f Z.„o„-rel. ^ ^ (5) 


where the tSZ spectral function (in intensity units) is 


2 (fcBTcMB)^ x^e^ 
(hc)^ {e^ — 1)^ 


X coth 




( 6 ) 


with X = hvflk^TcMB)- In the relativistic case, correc¬ 
tions to the spectral function arise that explicitly depend 
on the electron temperature: 

^ jtSZ.rel. [I + (ft, Te)] y (ft, M, z) . 

(7) 

We include relativistic corrections up to third order in 
k-BTe/irrieC^) following Ref. [34]. Note that relativistic 
corrections involving the cluster velocity also appear in 
general, but when averaged over all clusters in the uni¬ 
verse, these cancel to lowest order due to the differing 
sign of the line-of-sight velocity from cluster to cluster. 
The neglected, next-order contribution corresponds to 
(y) « 10“® [43] and its minor effect is left to future work. 
We also assume the contribution to (y) generated by the 
CMB dipole, (y)dipoie — 2.6 x 10“^, is subtracted using 
external measurements of the dipole amplitude [U]. 

We now briefly summarize the model used to calculate 
Eq. Q, with explicit details given in Ref. [16] and B15. 
We assume the lUAT4P5-|-eCMB-l-BAO-l-iTo maximum- 
likelihood cosmological parameters (e.g., cts = 0.817 and 
Vim = 0.282) [3^ and the halo mass function of Ref. [46] . 
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FIG. 1: The mean tSZ signal of the universe. The dashed 
cyan and solid blue curves show non-relativistic and relativis¬ 
tic calculations, respectively. The signal is dominated by hot, 
free electrons in galaxy groups and clusters. Error bars are 
shown only on the relativistic curve for clarity, and include the 
PIXIE instrumental noise, component separation noise, and 
cosmic variance (CV). PIXIE can detect the signal at 1470 ct 
significance; CV reduces the effective signal-to-noise to 230 (t. 
For comparison, the thin green curve shows the /r distortion 
(multiplied by 100 to render it visible). 


We use the ICM electron pressure profile fitting func¬ 
tion from Ref. m, which is extracted from cosmological 
hydrodynamics simulations [53]. This pressure profile 
matches a wide range of recent tSZ and X-ray obser¬ 
vations (e.g., [48U52] ) and fully specifies the ICM elec¬ 
tron pressure as a function of cluster mass, redshift, and 
cluster-centric distance. Finally, to calculate the rela¬ 
tivistic corrections, we use the Te(M, z) relation from 
Ref. [55], with a -1-20% correction applied to the masses 
derived from X-ray data in that work, to account for de¬ 
viations from hydrostatic equilibrium in the ICM |14j . 

The sky-averaged tSZ signal also receives contribu¬ 
tions from electrons in the ICM and during reioniza¬ 
tion. Our reionization model is described in Ref. |54] 
and B15. The IGM and reionization contributions are 
subdominant to the ICM signal by more than an or¬ 
der of magnitude. They are approximated well by the 
non-relativistic tSZ spectrum due to the electrons’ low 
temperature {k^Te few eV), and thus are fully char¬ 
acterized by the Compton-y parameter. We add these 
contributions to that from the ICM to obtain the total 
sky-averaged tSZ signal. We verify the accuracy of our 
analytic calculations by comparing to numerical simula¬ 
tions, hnding that the predicted (y) values agree to within 
2% (see BI5). 

Fig.g shows the mean tSZ signal of the universe, for 
both the non-relativistic and relativistic cases, as well as 
the /i distortion signal for comparison. We show below 



FIG. 2: Difference between relativistic and non-relativistic 
predictions for the mean tSZ signal of the universe. The solid 
magenta curve shows the difference between the tSZ predic¬ 
tions from Fig. [^ The dashed orange curve shows an ap¬ 
proximation based on moments of the optical depth-weighted 
IGM electron temperature distribution, which matches the 
full calculation to < 0.1% precision. The dash-dotted black 
curve shows the lowest-order “residual” (non-y/non-p) dis¬ 
tortion [^. The relativistic tSZ signal is partially degener¬ 
ate with this distortion, and thus makes it harder to access 
the r-type signal. The shaded blue area shows the PIXIE 
instrumental noise plus component separation noise, while 
the shaded red area shows the additional uncertainty from 
GV (note that important off-diagonal contributions are not 
shown). PIXIE can distinguish between the relativistic and 
non-relativistic predictions at 30a significance. 


that currently proposed experiments can measure the rel¬ 
ativistic effects at high precision. The non-relativistic re¬ 
sults can be summarized fully by the Compton-y param¬ 
eter. We find (y)icM = 1-58 x 10“®, (y)iGM = 8.9 x 10“®, 
and (y) reion = 9.8 X 10“® for the contributions from the 
ICM, IGM, and reionization, respectively. Note that 
{y)iCM depends sensitively on erg, the amplitude of mat¬ 
ter density perturbations (going roughly as dg); if we 
assumed Planck 2015 cosmological parameters [55] in¬ 
stead of WMAP9^ the prediction would be « 10% higher. 
Regardless, the ICM contribution dominates over those 
from the IGM and reionization, although it may be pos¬ 
sible to isolate the latter by masking the ICM using deep 
galaxy or cluster catalogs, or via cross-correlation tech¬ 
niques. All (y) contributions are much larger in ampli¬ 
tude than the y distortion signal. In agreement with 
early estimates [55], the total (y) is roughly one order of 
magnitude below the COBE-FIRAS bound. 

To emphasize the relativistic effects, Fig. [^ shows the 
difference between the non-relativistic and relativistic 
predictions. The fractional difference is « 1% &,i v < 500 
GHz, but is larger at higher frequencies (though the 
overall tSZ signal is smaller there). Fig. also shows 














































4 


the lowest-order “residual” (r-type) distortion (see Fig. 4 
of [56]) assuming an amplitude Hr = 10“®, which could 
be generated by (e.g.) decaying particles in the early 
universe (note that in our fiducial model Hr = 0; it is 
shown here only for reference).^ Fig. shows that the 
relativistic tSZ signal projects strongly onto the r-type 
distortion and thus renders its extraction more difficult. 

To interpret the relativistic tSZ signal, we consider a 
model based on moments of the optical depth-weighted 
ICM electron temperature distribution |35j : 


{{ksTeDr = 


[r) 


dz 


ICM 


d^V 

dzdn 


,,, dn 
dM-— X 

dM 


d^ntTT J dl Hein,!, M,z) [kBTe{n,l,M,z)]^ ,{8) 


where (t)icm is the mean Thomson optical depth in 
groups and clusters: 

JJdlne{di,l, Mys). ( 9 ) 

We make the approximation that only electrons in groups 
and clusters are hot enough to require a relativistic treat¬ 
ment, i.e., the IGM and reionization contributions are 
fully characterized by {y)iGM + (y) reion (which cannot be 
individually distinguished). We implement the moment 
formalism for the ICM electrons following Ref. |3S], con¬ 
sidering moments up to n = 4. 

The model thus contains six parameters: {{kBTe)'^)r 
for n = 1 to 4, (i/)igm + (2/)reion, and the total (y), 
which is interchangeable with (r)icM via (r)icM = 
iiy) - (y)lGM - (2/)reion)WeC^/(^Bre)r- The fiducial val- 
ues are ((/cBTe))^ = 0.208 keV, ((fcBre)^)^ = 0.299 
keV^, {{kBTe)^)r = 0.892 keV^, ((fcBTe)^)^ = 4.02 keV^ 
{y)iGM + (y)reion = 1-87 X 10“'^, and (y) = 1.77 x 10“® 
(equivalently, (t)icm = 3.89 x 10“^). The prediction of 
this model for the relativistic deviation is shown in Fig.[^ 
It agrees with the full calculation to < 0.1% precision 
over the entire frequency range. Thus, a measurement of 
the mean tSZ signal can be robustly interpreted in terms 
of these parameters. B15 presents a detailed breakdown 
of the mass and redshift contributions to these quantities. 

Prospects for Upcoming Experiments — PIXIE is a pro¬ 
posed Explorer-class mission designed to measure the 
absolute intensity and linear polarization of the sky at 
1.6 degree resolution in 400 linearly-spaced frequency 
channels from 30 GHz to 6 THz [T5j. PIXIE's un¬ 
precedented spectral coverage and resolution will al¬ 
low for exquisite component separation, yielding (low- 
angular resolution) maps of the tSZ signal. Galactic dust 


^ Note that the r-type spectral shape was derived while neglecting 
relativistic tSZ effects (i.e., considering only pure y and p) [56]. 


and synchrotron emission, and the cosmic infrared back¬ 
ground, as well as precise constraints on primordial ten¬ 
sor fluctuations from the inflationary epoch. We con¬ 
sider 80 frequency channels of width Xvp = 15 GHz 
ranging from 30 GHz to 1230 GHz, with top-hat band- 
passes. The tSZ signal is negligible at higher frequen¬ 
cies, but we assume the PIXIE data in these channels 
will be used to measure dust foregrounds with sufficient 
precision to subtract them to the level of the instru¬ 
mental noise at the lower frequencies (see Sec. 3.2 of 
Ref. [H]). The instrumental noise in each channel is 
A/p = 5 X 10“^® Wm“^ Hz“^ sr“^, which we take to 
be uncorrelated between channels [El [56]. To account 
for noise arising from component separation, we model 
the sky-averaged signal using the terms in Eq. (!]). In 
the Fisher matrix approximation, we find that marginal¬ 
izing over the AT and /i components yields an effec¬ 
tive total noise on the tSZ contribution of A/p^g^ = 
8.8 X 10“^® Wm“^ Hz“^ sr“^, implying a raw detection 
significance of 1470cr for the sky-averaged tSZ signal 
shown in Fig. [^ 

In addition to the diagonal contribution to the channel- 
channel covariance matrix from instrumental noise, there 
is a (non-diagonal) contribution from cosmic variance: 


Gov 


i/v' _ 

cv — 


C w' 

47r/sky ’ 


( 10 ) 


where is the zero mode of the relativistic tSZ cross¬ 
power spectrum at frequencies v and v', and /sky is the 
observed sky fraction. We derive this result in B15. We 
use the exact, full-sky results for the tSZ power spec¬ 
trum from Ref. m to compute accounting for rel¬ 

ativistic effects using the same model as in our mean tSZ 
calculation. We assume /sky = 0.75 for PIXIE, which 
is reasonable given its spectral coverage and ability to 
remove foregrounds. The full covariance matrix is 

Cov'^^'' = + Coy'SC , (11) 

where 6'^'' is the Kronecker d-function. The CV term 
contains non-negligible off-diagonal contributions, since 
nearly the same population of groups and clusters is re¬ 
sponsible for the signal at different frequencies. 

Including the CV contribution to the covariance ma¬ 
trix in Eq. 0, the effective significance of the sky- 
averaged tSZ measurement decreases to 230cr, which is 
still extremely promising. The ultimate CV limit for a 
full-sky, zero-noise experiment is ~ lOOOcr, depending on 
how much information in the high-frequency relativistic 
tail can be accessed. It would be possible to decrease the 
CV by masking massive, low-redshift clusters that con¬ 
tribute significantly to but not as significantly to 

the mean tSZ signal, similar to low-/ measurements of 
the tSZ power spectrum mm- 

Fig. [2 shows error bars computed from the square root 
of the diagonal elements of the covariance matrix, includ- 
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ing CV. The error bars are difficult to see, but are more 
visible in Fig.[^ where we show the effective PIXIE noise 
(including component separation degradation) and the 
full errors including CV (note that the off-diagonal CV 
contributions are important, but not shown on the plot). 
Using the full covariance matrix, we find that PIXIE can 
distinguish between the relativistic and non-relativistic 
mean tSZ predictions at 30 ct significance. 

We use the Fisher matrix formalism to forecast PIXIE 
constraints on the parameters of the moment-based 
model presented in Fig. Degeneracies prevent mean¬ 
ingful constraints from being placed on most of the pa¬ 
rameters at pixie’s noise level, but if only (y) and 
{kBTejr are allowed to vary, we find marginalized Icr un¬ 
certainties of = 8.3 X 10“^ and cr{kBT^)^ = 7.6 x 10“^ 
keV. These values correspond to fractional uncertainties 
of 0.5% and 4%, respectively. PIXIE will thus improve 
the existing COBE-FIRAS bound on (y) by three or¬ 
ders of magnitude, while simultaneously measuring the 
mean optical depth-weighted ICM electron temperature 
at high precision. It may be possible to further improve 
these forecasts by combining with constraints on (r)icM 
from patchy r estimators applied to CMB polarization 
maps [5^. Finally, note that the inclusion of (fcBTe)r 
degrades the constraint on /r by « 10% compared to a 
forecast considering only the non-relativistic tSZ signal 
(/i has been marginalized over in the constraints quoted 
above). 

Outlook — Spectral distortions of the CMB contain a 
vast and complementary set of cosmological information 
to that contained in the temperature fluctuation power 
spectrum and other spatial statistics. Amongst this in¬ 
formation lies the mean tSZ signal of the universe, which 
can be detected at a raw significance of nearly 1500cr with 
PIXIE, yielding a precise constraint on the total thermal 
energy in electrons in the observable universe. 

In B15, we translate our forecasts into direct con¬ 
straints on parameters describing physical models of the 
ionized gas. In addition, the relativistic effects considered 
here lead to mixing between the low-redshift tSZ signal 
and the primordial y and “residual” (r-type) spectral dis¬ 
tortions, which must be accounted for in future analyses 
(see Fig.[^. Overall, a measurement of the mean tSZ sig¬ 
nal of the universe represents an important step forward 
in our understanding of the thermodynamics of ionized 
gas in the universe. 
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